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Abstract 

When a two-level quantum dot and a plasmonic metal nanoantenna are resonantly coupled by the 
electromagnetic near field, the system can exhibit a Fano resonance, resulting in a transparency dip 
in the optical spectrum of the coupled system. We calculate the nonlinear response of such a system, 
for illumination both by continuous-wave and ultrafast pulsed lasers, using both a cavity quantum 

o 

O electrodynamics and a semiclassical coupled-oscillator model. For the experimentally relevant case 

,__! of meV thermal broadening of the quantum-dot transition, we predict that femtosecond pulsed 

> 

illumination can lead to a reversal of the Fano resonance, with the induced transparency changing 

m 

I/*} into a superscattering spike in the spectrum. This ultrafast reversal is due to a transient change in 

—f the phase relationship between the dipoles of the plasmon and the quantum dot. It thus represents 

O 

a new approach to dynamically control the collective optical properties and coherence of coupled 






nanoparticle systems. 
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A hybrid system of a semiconductor quantum dot (QD) in the near field of a plasmonic 
metal nanostructure can exhibit qualitatively different optical properties than its individual 
components [IH3|. For sufficiently strong resonant coupling between the QD exciton and 
the plasmon, the optical spectrum can exhibit Fano interference jH [5]: the QD creates a 
dramatic "dipole-induced transparency," [6] suppressing absorption and scattering in spite 
of its relatively meager oscillator strength. A classical model can describe this effect in 
the linear response limit. Calculating the nonlinear response of the hybrid system requires 
that the QD, at least, be modeled quantum mechanically. Semiclassical (SC) models that 
treat the QD as a two-level system have predicted novel nonlinear-optical effects, including 
a "nonlinear Fano effect" and optical bistability [U ITHTO]. Treating both plasmon and 
QD quantum mechanically further refines the picture [TTHT5] . predicting in particular a 
suppression of induced transparency, and of bistability, due to additional dephasing not 
accounted for in the SC model [131 HB]- 

These phenomena have all been predicted in the regime of very narrow QD linewidths, 
on the order of 10 /xeV, corresponding to liquid-helium temperatures. Although this can be 
realistic when QDs are coupled to practically lossless components such as photonic-crystal 
cavities, absorptive heating may render such temperatures infeasible for QDs coupled to 
plasmonic nanostructures. We therefore consider here the more experimentally achievable 
parameter regime of meV QD linewidths. 

In this regime, we predict a new nonlinear phenomenon for femtosecond pulsed excita- 
tion; namely, for particular fluences the Fano resonance reverses, resulting in a coherently 
enhanced cross section rather than an induced transparency. This is a dynamical analog of 
the transition from electromagnetically induced transparency to superscattering |ITj , arising 
from a transient change in the phase relationship between the QD and the plasmon. It is 
thus fundamentally different from the nonlinear Fano effect, a steady-state change in the 
Fano lineshape arising from the dependence of the QD-plasmon coupling strength on the 
incident field intensity [U [12] . 

Our treatment starts with a quantum-mechanical model of the hybrid system. We fol- 
low a previously developed cavity quantum electrodynamics (CQED) formalism (H [Hi [13], 
extending the previous work to the regime of broader QD linewidths and to consideration 
of the transient optical response. The underlying basis states are \qs), where q G {0, 1} 
indexes the QD energy levels and s € { 0, 1, 2, . . . } indexes plasmon energy levels. Lowering 



and raising operator pairs for the QD and plasmon are (a,a + ) and (b,b + ), respectively. 
The dipole operators are then fi q = d q (a + a + ) and p s = d s (b + b + ), where d q and d s are 
the transition dipole moments of the QD and plasmon, respectively, and the total dipole 
operator is p = p s + [2 q . The evolution of the density operator p(t) is governed by 

where H is the system Hamiltonian, and L(p) is a Lindblad superoperator accounting for 
disspiation. Specifically, H = Hg+Hs+Hi+Hj, where H q = hu q a + a is the QD Hamiltonian, 
H s = hw s b + b is the plasmon Hamiltonian, Hi = —hg(a + b + ab + ) describes plasmon-QD 
coupling, and Ha = —E(t)p describes driving by an incident electric field E{t) [18]. The 
Lindblad superoperator is [TT] 

ru. „ „ * A 

L q (p) = — — {a^ap + pa* a — 2d per') — 72 (a^ap + pa^a — 2a^apa^a) . (2) 

The first term describes spontaneous emission with rate 71 = T£" , and the second term 
describes dephasing with rate 72 = T 2 _1 . In order to solve Eq. (Ill), we define a maximum 
number of plasmons, N s , above which the density matrix elements are negligible. Solution 
then involves integrating 0(N%) coupled ordinary differential equations, or, at steady state, 
O(N^) coupled algebraic equations [IB]. Once a solution is obtained, the total dipole is 
calculated according to p(t) = Tr (p(t)p). 

A computationally simpler approach is a SC or Maxwell-Bloch model, in which the plas- 
mon dipole, p s (t), is treated classically and the QD is treated with Bloch equations [19] or 
their generalization [20J for its reduced density matrix, p® D (t): 
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(i s + 7 s /i s + uj 2 s p s = A S [E + Jp q ] , 



Pi = W g p 2 - 72Pl, P2 = -UJqPl ~^-[E + Jp s ] p 3 - 72P2 , (3) 

9d 

p 3 = -^iE + Jp s }p 2 - 1 s 1 c (p 3 + l), 

where E = E(t), pi = 2Rep^ D (t), p 2 = -2Imp^ D (t), p 3 = p^ D (t) - Po D (t), and p q = d q p 1 . 
One can relate A s to the CQED parameters by solving for the steady state of p s with J = 
[H], giving A s = Au s d 2 s /h. Comparing classical and quantum dipole interaction energies 
gives J = hg/(d q d s ) [H]. Eqs. [3] can be numerically integrated with or without the rotating 
wave approximation [H]; the results are essentially identical for the present problem. The 
SC model involves the solution of a fixed, small number of ordinary differential equations 



regardless of intensity, unlike the CQED model, which becomes computationally costly for 
large N s . 

A flaw in the SC model is that an excited QD in the dark cannot couple to the plasmon, 
but does so in CQED via the Purcell effect (TTJ HE] 121] . With u q = u s and 7 S ^> 71,72, 
and g, the QD's effective decay constant due to the Purcell effect is [18] jf 5 ~ 4g 2 /j s . We 
therefore term the SC model with 7^ = 71 as the "naive" SC model, and with 7^° = 7^ 
as the "corrected" SC model. 

Having obtained //(£) using either the CQED or SC model, the absorption cross section 

is found according to jTHl E2] : 

k 
cr a bs(u) = — Im[a(u)] , (4) 

where k = s/e^dcu/c, with e me d being the relative dielectric constant of the surounding 
medium, and 

«M = — — r — : -L, ^ \ > (5) 

where integration is over an optical cycle after steady state is reached or over all time for 
pulsed excitation. While we discuss absorption cross sections here, scattering spectra exhibit 
nearly identical trends. 

The model parameters are obtained by fitting to spectra for a realistic system calculated 
with the discrete dipole approximation (DDA) [181 E3]. As illustrated in Fig. [II two Au 
prolate spheroids with semi-major and semi-minor axes of 15 nm and 10 nm, respectively, 
are arranged coaxially with a gap of 6 nm. A 4 nm diameter CdSe QD is placed in the 
center of the gap. The system is embedded in a medium with dielectric constant e me d = 2.25, 
typical of a polymer or glass. The QD dielectric constant is taken to be a Lorentzian function 
with center frequency chosen to match the plasmon frequency of the metal nanostructure, 
and with linewidth corresponding to temperatures of 50 - 100 K [5] [HI El]. The fitting 
gives hu s = 2.042 eV, d 8 = 2990 D, h-? s = 150 meV; hu q = 2.042 eV, d q = 13.9 D, 
/ry 2 = 1.27 meV; and hg = 10.8 meV [18]. As seen in Fig. [IJ the CQED and SC results 
are in excellent agreement with each other and in good agreement with the DDA spectrum. 
The QD spontaneous emission rate is calculated according to 71 = u^y/e me dd^/ (3iie hc 3 ) 
[25], giving foyi = 268 neV or T\ = 7^ x = 2.46 ns; this, in turn, gives h^ = 3.02 meV, 
corresponding to T^ =218 fs. 

We begin by considering steady-state spectra. Fig. [2] shows that the Fano resonance dip 
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FIG. 1. Linear absorption spectrum of an Au-CdSe-Au hybrid nanoparticle system (illus- 
trated in the inset), calculated using the discrete dipole approximation (DDA), cavity-quantum- 
electrodynamics model (CQED), and semiclassical model (SC). 

disappears as the incident intensity is increased, due to saturation of the QD transition. The 
corrected SC and CQED models are in excellent quantitative agreement for high and low 
applied fields. The naive SC model is gravely in error, but can be brought into agreement 
with the corrected SC model by multiplying the incident intensity by the Purcell factor 
[TU [15]. This validates the use of the SC model, which is particularly important for simula- 
tions at the high intensities for which CQED calculations are computationally prohibitive. 
The results for the current system contrast with predictions for systems with narrow QD 
linewidths; for these systems, the SC formulation gives a deeper transparency even in the 
linear regime, due to quantum-optical dephasing that is ignored by the SC model |H]. In 
our system, thermal dephasing dominates over this vacuum-field-induced dephasing. As the 
field is increased, however, the quantum-optical dephasing increases and eventually becomes 
comparable to the thermal linewidth, resulting in a small disagreement between our CQED 
and SC predictions at moderate fields. 

We next consider the system response to a Gaussian pulse, E(t), that has a 20 fs full 
width at half the maximum intensity and a center frequency, hu> = 2.042 eV, resonant 
with the plasmon and exciton. Transient spectra obtained by Fourier transformation of the 
time-domain response are shown in Figure [3j At low intensities, the transient spectrum is 
identical to the linear steady-state spectrum. Strikingly, however, at certain higher ffuences 
the resonant dip reverses to form a narrow spike. The corrected SC model remains in 
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FIG. 2. (a) Steady-state absorption spectra calculated using the CQED (dots) and corrected 
SC (solid) model for various incident intensities. Successive spectra are displaced vertically by 
2 x 1CT 11 cm 2 , (b) Intensity dependence of steady-state absorption cross-sections calculated using 
the CQED, naive SC and corrected SC steady-state models, at a photon energy of 2.042 eV. 



excellent agreement with the CQED model. As the pulse fluence increases beyond the range 
that is readily computable using the CQED model, the SC model predicts recurrences of 
the ultrafast reversal suggestive of Rabi oscillations. 

The ultrafast reversals of the Fano resonance arise from the transient phase relationships 
of the plasmon and QD dipoles with respect to the incident light, as illustrated in Fig. |4} At 
the beginning of the incident pulse, as at steady state, the plasmon lags the driving laser by 
the 7r/2 phase difference expected for a resonant oscillator [15]. The QD is driven primarily 
by the plasmon and thus lags the laser by an additional tt/2 phase difference. At low fluence, 
this phase relationship continues until the pulse is complete and the short-lived plasmon has 
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FIG. 3. (a) Absorption spectra calculated using the CQED (dots) and corrected SC (solid) models, 
for ultrafast pulsed excitation. Successive spectra, corresponding to varying fluences, are displaced 
vertically by 2 x 10~ n cm 2 . (b) Fluence dependence of the absorption cross-section for ultrafast 
pulsed excitation, calculated using the CQED and SC models, at a photon energy 2.042 eV. 



decayed. Then, the longer-lived QD, still oscillating with a ir phase lag relative to the laser, 
drives the plasmon; the plasmon thus acquires a it + (tt/2) = (37r/2) phase lag, partially 
canceling its earlier oscillations in the spectral domain and producing the observed linear 
Fano dip. By constrast, at higher fluences, the QD population exhibits Rabi oscillations, 
reaching unity and then coherently being driven back down. This reverses the sign of the 
QD dipole [2H], so that the lag of the QD phase relative to the laser is now zero. The phase 
frustration that previously led to transparency is replaced by a constructive interference 
that leads to induced absorption or superscattering. 

Ultrafast reversal is thus due to a change in the phase of the coherent interaction between 
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FIG. 4. Time evolution of the QD-plasmon system under pulsed excitation, calculated using the 
CQED model for two fluences. (a,d) Pulse's electric field and QD population. (b,e) QD's dipole 
and phase lag relative to the pulse. (c,f) Plasmon's dipole and phase lag. Dashed lines indicate 
times at which a phase jump occurs. 



the QD and the plasmonic metal nanostructure. A similar change has been demonstrated 
in plasmonic or metamaterial systems that exhibit Fano resonances due to the coupling 
between bright modes and dark modes [27] . In this case, the sign of the interference can 
be controlled through careful selection of the linewidths and coupling strengths [28], or by 
adding a retardation-based phase delay [29] . Similarly, a change from steady-state trans- 
parency to enhanced absorption has been predicted in coupled QD-plasmon systems by 
changing the size of the metal nanoparticle and the detuning [15j. In these systems, reversal 
from destructive interference to constructive interference can be controlled only statically, by 
changing the structure of the system. In our QD-plasmon system, by contrast, the reversal 
can be controlled dynamically, by changing the fluence of the incident ultrafast pulses. 

Not all QD-plasmon systems that exhibit Fano resonances will exhibit ultrafast reversal. 
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The reversal requires that the QD dipole oscillates significantly longer than the plasmon's 
intrinsic lifetime. In other words, we require 7^ ff -C j s , which implies g <C 7 s /2. However, 
g must also be large enough to give a Fano resonance, which requires g > ^/j s j2 /4 [TB] . 
These constraints are more easily satisfied for 72 ^C 7s- The large 7 S afforded by a lossy 
plasmonic component points to the intrinsically plasmonic nature of ultrafast reversal: such 
constraints are unlikely to be satisfied by a high-finesse photonic resonator. 

Sufficiently narrow QD linewidths can be obtained by cooling to liquid-nitrogen temper- 
atures, even if we account for absorption-induced heating of the plasmonic nanostructure. 
Achieving the required coupling strengths is a greater experimental challenge, but should be 
feasible using chemically synthesized components and directed assembly, such as DNA-based 
assembly of colloidal QDs and metal nanoparticles [30] . The availability of low-cost assem- 
bly methods is an advantage for these systems as compared to traditional CQED systems, 
which require complex and expensive top-down fabrication 



Our treatment of the optical response of coupled QD-metal nanostructure systems has 
employed a larger thermal dephasing rate for the QD than has generally been considered 
in previous treatments. Although this means that certain phenomena requiring a high 
degree of coherence are suppressed, significant quantum-optical effects remain. First, the 
saturation of the Fano resonance is the principal optical nonlinearity at steady state, and the 
intensity at which this saturation occurs is due to a balance between two inextricable aspects 
of QD-plasmon coupling: plasmonic field enhancement lowers the incident fields required 
for saturation, while the Purcell effect increases the required fields. Second, we predict 
that the Fano resonance can undergo a reversal, changing from a transparency dip into a 
superscattering spike, when excited by femtosecond laser pulses with appropriate fluence. 
This ultrafast reversal represents a new means to coherently control optical interactions 
among nanostructures. 
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